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ABSTRACT 


A  mixed  layer  dynamic  model  fcr  the  structure  and  evolution 
cf  aercscls  in  a  marine  regime  is  presented.  The  aerosol 
spectrum  is  divided  into  continental  and  marine  components, 
with  a  constant  vertical  continental  profile  assumed;  bcth 
in  and  above  the  mixed  layer.  The  aerosol  spectrum  is  trans¬ 
formed  into  a  reference  relative  humidity  of  8091.  The 
temporal  evolution  of  the  aerosol  spectrum  is  predicted  from 
rate  equations  which  require  a  specification  of  the  surface 
production  rate,  the  entrainment  rate  (We),  and  the  mixed 
layer  depth  (h) .  The  model  was  tested  against  the  data  set 
obtained  in  the  Honterey  Bay  during  the  HAGAT  80  experiment. 
The  model  was  initialized  with  both  observed  HAGAT  data  and 
an  equilibrium  initial  value,  gathered  from  the  JASIN 
Experiment.  The  model  was  run  for  radii  equal  to  0.8,  2.0, 
5.0  and  10.0  microns.  The  significance  between  the  observed 
initial  values  and  the  JASIN  data  is  that  the  HAGAT  data 
were  observed  from  the  same  air  mass  as  the  initial  atmos¬ 
pheric  data  and  verification  aerosol  data.  The  JASIN  initial 
aerosol  data  are  based  on  an  equilibrium  state  as  a  function 
of  only  wind  speed  and  reference  relative  humidity.  The 
model  continuously  generated  a  correct  gain  or  loss  of 
aercscl  concentrations  as  defined  by  the  observed  HAGAT 
data,  and  in  most  cases  the  model  output  is  within  one  crder 
cf  magnitude  of  the  observed  values. 
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I.  IH1BCD0CTI0H 


There  is  an  increased  emphasis  in  the  ability  to  fore¬ 
cast  the  behavior  and  evolution  of  aerosols.  The  scattering 
and  absorption  of  light  are  influenced  by  aerosol  distribu¬ 
tions  in  the  atmospheric  boundary  layer  and  affect  the 
performance  of  optically  guided  weapon  systems.  The  Air 
Force  is  particularly  interested  in  how  aerosol  extinction 
affects  the  use  of  precision  guided  munitions  (PGM) 
/Cottrell  et  al, 1979 ]  .  The  Department  of  Defense  (COD)  has 
PGM  *  s  that  operate  at  differing  wavelengths  which  range  from 
the  visible  to  the  microwave  regions.  PGM's  have  a  greater 
ability  to  hit  a  target  than  conventional  munitions; 
however,  the  controlling  factor  is  the  ability  of  the 
guidance  system  to  "see"  the  target.  This  ability  is  depen¬ 
dent  on  the  wavelength  for  which  the  sensors  are  designed 
and  the  properties  in  the  intervening  atmosphere.  The 
degrading  properties  in  the  atmosphere  are  principally 
molecular  absorption  and  aerosol  scattering.  The  wavelengths 
for  the  guidance  systems  are  designed  such  that  the  molec¬ 
ular  absorption  is  minimized:  therefore,  scattering  by  aero¬ 
sols  becomes  the  main  concern  once  a  suitable 
absorption- free  window  has  been  selected. 

The  model  under  consideration  includes  the  behavior  of 
marine  aerosol  components  as  well  as  previous  continental 
components.  Estimating  the  influences  of  marine  aerosols  on 
elect zo-cptical  (EO)  systems  has  been  studied  /Barnhart  and 
Streets,  1970/  .  Particle  sizes  of  interest  are  those  asso¬ 
ciated  with  locally  generated  sea  salts  because  of  their 
effects  on  IR  as  well  as  visible  wavelengths.  The  size 
distributions  of  sea  salt  particles  show  a  variance  of 
several  orders  of  magnitude. 
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The  ability  to  forecast  the  behavior  of  aerosol  extinc¬ 
tion  fron  synoptic  scale  patterns  would  help  in  the  deciding 
|  '  which  type  cf  weapon  system  to  employ.  Because  some  systems 

are  launched  from  the  air,  it  is  important  that  such  a 
profile  include  the  vertical  distributions  of  aerosol  parti¬ 
cles.  Models  exist  for  estimating  vertical  extinction 
|  profiles,  but  they  have  not  been  sufficiently  verified.  To 

do  this,  profiles  cf  actual  aerosol  data  must  be  gathered 
and  compared  with  the  model  forecasts. 

Mcdels  in  current  use  are  based  on  parameterizaticns  of 
the  effects  of  relative  humidity  and  wind  speed  on  the  equi¬ 
librium  aercsol  distributions  /wells,  et  al,  1977^  .  Recent 
evaluations  have  shown  that  these  models  are  limited  tc  mean 
distributions  (i.e.,  the  average  aerosol  concentration  at  a 
given  wind  speed  and  humidity)  /Fairall  et  al,  1982a/  . 
Models  are  limited  because  seme  processes  in  the  atmospheric 
mixed  layer  which  affect  aerosol  concentrations  are  not 
considered,  namely  entrainment  and  inversion  height  changes. 

The  purpese  of  this  study  is  to  present  and  evaluate  a 
model  which  includes  the  meteorological  processes  that 
adequately  describe  the  whole  marine  atmospheric  boundary 
layer  (NABL) .  in  inversion  represents  a  cap  to  the  vertical 
transport  of  surface  generated  aerosols,  and  is  not 
accounted  fer  in  previous  models.  The  top  of  the  boundary 
layer  is  capped  by  the  marine  inversion,  where  entrainment 
of  overlying  air  takes  place.  Because  entrainment  mixes 
clear  (nen-marine)  air  into  the  marine  layer,  this  process 
is  as  important  as  surface  layer  fluxes  in  determining  equi- 
libriui  concentration.  This  entrainment  process  is  included 
in  the  mcdel. 

Evaluaticn  of  the  icdel  cutput  will  be  done  with  the  data 
from  at  experiment  entitled  Marine  Aerosol  Generation  and 
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Transput  (MAGAT)  .  The  experiment  was  conducted  in  the 
vicinity  of  Monterey  Eay,  Car  during  the  period  cf  28  April 
to  9  Kay  1980.  The  purpose  of  the  experiment  was  to  examine 
the  compatibility  of  optical  and  micrometeorological • propa¬ 
gation  theory,  and  tc  extend  dynamic  models  of  the  evolving 
MABL  tc  include  aerosol  and  turbulence  profiles 
/fairall,  1980  and  Fairall  et  al,1980/  .  Two  platforms,  the 
E/V  ACANIA,  and  an  aircraft  were  used. 

In  this  study,  twc  24-hour  periods  were  chosen  for  evalu¬ 
ation;  3  Bay,  beginning  at  1200  PDT  and  6  May,  beginning  at 
1800  PET •  These  periods  were  chosen  because  the  boundary 
layer  was  undisturbed  (no  fronts  closer  than  100  nm)  for  24 
hours  prior  to  the  starting  times.  Additionally,  the  times 
were  selected  dee  tc  the  proximity  of  the  aircraft  and 
surface  ship  during  the  experiment  times.  The  approach  was 
to  describe  the  synoptic  conditions  from  24  hours  in  advance 
cf  the  model  forecast  through  the  end  of  each  forecast,  and 
compare  the  model  output  with  the  actual  MAGAT  findings. 


II.  BiCKGHOOHD 

1.  SESCB1PTION  OF  EXISTING  BODELS 

Current  models  for  estimating  aerosol  equilibrium  distri¬ 
butions  use  meteorological  inputs  of  (10-m  height)  wind 
speed  atd  relative  humidity.  These  two  quantities  are 
considered  because  of  their  role  in  generation  and  transport 
(wind)  and  growth  (humidity)  of  aerosols.  The  Navy's 
lells-Munn-Katz  (WMK)  /wells  et  al,  1977/  is  an  example  of 
this  concept.  The  performance  of  this  type  of  model  has 
teen  studied  with  data  obtained  in  the  northern  Atlantic  and 
eastern  Pacific  Ocean  areas.  The  model  output  compared  with 
these  data  is  shewn  in  figure  2.1,  depicting  a  height  depen¬ 
dence  cf  tctal  aeroscl  volume  from  a  sample  set  of  eastern 
Pacific  data.  These  profiles  correspond  to  (1)  the  observed 
sea  salt  volume,  V;  (2)  the  observed  sea  salt  volume 
adjusted  to  80*  relative  humidity,  V0  ;  and  (3),  the  HHK 
predicted  volume  adjusted  to  80*  relative  humidity 
(circles),  it  is  clear  from  this  figure  that  within  the 
mixed  layer  the  observed  decrease  of  aerosol  volume  with 
height  is  less  than  the  model  predicts.  The  surface  gener¬ 
ated  aerosols  appear  to  be  well  mixed  below  the  inversion 
when  normalized  to  80%  relative  humidity. 

The  assertion  that  existing  models  can  predict  only  a 
mean  value  appears  in  the  results  obtained  in  the  North 
Atlantic,  from  the  JASIN  Experiment.  Figure  2.2  compares  a 
single  radius  size  cf  five  microns  and  the  corresponding 
model  prediction.  The  values  and  trends  in  the  predicted  and 
mean  results  are  in  reasonable  agreement.  However,  the 
standard  deviation  is  three  times  the  mean;  if  one  assumes  a 
normal  distribution,  only  67%  cf  the  observed  aerosol 
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figure  2.1  Height  dependence  of  aerosol  voluae  and  relative 

naiidlty  above  the  ocean. 
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distributions  will  be  wit  fair 


factor  of  three  of  the 


average.  This  comparison  emphasizes  the  point  that  nc  matter 
how  accurately  a  model  predicts  an  average  aeroscl  density 
at  a  given  wind  speed  and  relative  humidity,  the  factcr  of 
three  BBS  variations  cannot  be  eliminated  without  consid¬ 
ering  mere  meteorological  parameters. 

B.  IBS  BBS  BOUHDABT  IATEH  BODE! 

The  distinguishing  feature  of  this  model  is  the  charac- 
terizaticn  of  the  HAEL,  which  is  convectively  mixed  up  tc  a 
height  h,  and  capped  by  an  inversion.  The  atmospheric 
profile  is  depicted  in  figure  2.3a,  representing  a  cloud- 
free  mixed  layer  where  water  vapor  mixing  ratio  (g)  ,  and 
virtual  potential  temperature  (ev)  are  "well-mixed",  ie, 
independent  of  height  in  the  mixed  layer.  The  assumed 
vertical  aerosol  profile  is  shown  in  figure  2.3b.  The  model 
produces  a  24-hcur  time  evolution  of  an  aerosol  spectrum, 
requiring  a  prediction  of  the  following  at  each  time  step: 

1)  surface  production  rate  of  marine  aerosols 

2)  entrainment  rate  at  the  top  of  the  mixed  layer 

3)  mixed  layer  depth. 

retails  cf  how  these  parameters  are  input  into  the  model  are 
discussed  in  chapter  IV.  Ike  mathematical  relationship 
between  the  time  rate  of  change  of  the  aerosol  volume  spec¬ 
trum,  dv/dr,  and  these  three  parameters  is  shown  in  the 
following  equation 


dv/dt  =  (<w’v'>0  -(We+Wkm)  v)/ h 
1  2  3 


(2.1) 


The  model  predicts  the  evolution  of  aerosol  at  five  radii 
(0.8,  2.0,  5.0,  10.0,  and  15.0  microns)  of  both  the  conti¬ 
nental  and  the  marine  (sea  salt)  components  in  the  mixed 
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.3  Profile  of  a)  Hired  layer 
vertical  aerosol  distribution. 


layer.  The  actual  time  evolution  of  each  radius  is  output  in 
the  fcra  of  dv/dr,  with  the  units  of  urn2  /cm3. 

The  model  assumes  rc  appreciable  concentration  of  marine 
aerosols  above  the  mixed  layer.  The  entrainment  of  air  into 
the  mixed  layer  will  not  cause  an  element  of  these  marine 
particles  tc  "escape"  into  the  free  troposphere.  The  ocean 
acts  as  a  source  for  marine  aerosols,  primarily  through  the 
generation  of  white  caps.  This  is  the  only  input  for  the 
marine  component.  The  ocean  also  acts  as  a  sink  for  the 
continental  aerosols,  which  are  generally  of  a  smaller 
radius.  The  entrainment  process  at  the  top  of  the  mixed 
layer  mixes  clear  (ncn-marine)  air  into  the  N&BL.  This 
process  could  be  as  important  as  the  surface  layer  fluxes  in 
deter lining  an  aeroscl  concentration. 

C.  BOOH  INPUTS  AND  AEBOSOl  INITIALIZATION 

The  model  is  designed  so  that  with  the  exception  of 
initializing  the  aerosols,  all  of  the  calculated  parameters 
are  tasec  cn  inputs  from  surface  based  cbservaticns,  and 
representative  soundings.  The  input  consists  of: 

1)  Non-meteorological  inputs  of  latitude,  julian  day,  and 
the  lccal  start  time.  These  are  used  for  a  diurnal  radia¬ 
tive  heating  /ceding  package.  Diabatic  warming  has  an 
obvious  impact  on  the  life  and  strength  of  the  mixed 
layer.  Additionally,  sea-surface  temperature  (SST)  is 
input  at  the  start  time.  Up  to  ten  forecast  SST  values 
can  be  input  during  the  24-hcur  period  as  well. 

2)  Surface  wind  speed  and  direction  at  the  start  time  and, 
as  with  SST,  up  to  ten  additional  forecast  values. 

3)  Pica  figure  2.3,  temperature  and  relative  humidity  are 
parameterized  by  virtual  potential  temperature  and  water 


vapcr  Tiling  ratio,  respectively.  Mixing  ratios  are  input 
above  the  nixed  layer,  as  well  as  the  lapse  rate  above  the 
"  jaxp” . 

An  initialization  of  the  six  radii  concentrations  in 
three  regimes  is  reguired.  To  do  this  it  is  necessary  to 
understand  the  actual  distributions  of  the  continental  and 
narine  ccnpcnents  in  the  nixed  layer.  For  the  purposes  of 
evaluating  the  nodel,  assunptions  must  be  made  cn  the 
distributiens  of  the  acquired  MAGAT  data.  The  aircraft 
instruments  could  net  distinguish  the  chemical  make-up  of 
the  irdividual  particulates;  therefore,  it  is  not  knewn,  for 
a  given  radii  concentration  in  the  mixed  layer,  hew  much  is 
continental  and  how  much  is  narine.  The  concentrations  of 
each  size  were  calculated  in  the  nixed  layer,  and  then  above 
the  nixed  layer.  Based  on  the  assumption  that  there  are  no 
xarine  concentrations  above  the  mixed  layer,  we  conclude 
that  the  concentrations  calculated  from  above  are  all  conti¬ 
nental,  and  therefore  linerally  subtracted  from  the  concen¬ 
trations  in  the  tixed  layer,  leaving  only  marine 
concentrations. 

A  further  adjustxent  of  the  input  aerosol  values  is 
required,  based  on  changes  in  relative  humidity  during  the 
data  collection  tine.  He  have  stated  that  the  growth  of 
aeroscls  is  a  function  of  relative  humidity.  Therefore,  a 
distribution  of  aeroscls  gathered  at  9051  humidity  cannct  be 
directly  compared  to  another  distribution  gathered  at  80  % 
humidity;  a  reference  humidity  is  required.  Consider  the 
aerosol  volume  spectrum 


v(r)  =  4/3«r3  n(r) 


(2.2) 
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where  n  (r)  is  the  natter  density  spectrum.  V  (r) 
as  the  volume  of  aerosol  particles  per  cm3  at 
saturation  ratio  of  S*80%.  i  humidity  growth 
jl airall,  et  al,  1982aJ  is  defined  as 


is  defined 
a  reference 
factor  G  (S) 


G(S)  =  .81  exp(. 0665/(1.58  -S)> 


There  are  further  considerations  concerning  advecticn. 
Clearly  all  aerosols  do  not  originate  locally;  both  above 
and  telow  the  inversion  aerosols  are  advected  into  a  local 
region.  In  terms  of  aerosol  density,  entrainment  acts  as  an 
aerosol  flux  out  of  the  boundary  layer  because  the  concen¬ 
trations  above  and  below  the  layer  are  different.  In  the 
model,  the  entrainment  acts  on  the  "jump"  across  the  inver¬ 
sion.  For  the  purpose  of  evaluating  the  the  evolution  of  a 
local  concentration  cf  aerosols,  the  model  will  neglect 
horimcntal  advection  and  further  assume  a  negligible  local 
production  of  the  ccntinental  component.  The  model  takes 
into  consideration  a  Stokes  gravitational  fallout  term, 
jsu,  1S7j§  .  The  fallcut  rates  above  and  below  the  inversion 

are  different  because  of  the  change  in  aerosol  spectra, 
caused  by  the  humidity  growth  factor.  The  Stokes  velocity 
^airall,  et  al,  1S€2a/  is  calculated  from  the  following 
equation  where  pw  is  the  density  of  the  droplet  and  £  is  the 
kinematic  viscosity  cf  the  air. 


Wk  =  2g(pw  -p)  r2  G2(S)  /  <9Ep) 


(2.4) 


£.  ICOIIIBFIOH  AS  BOS  CL  BOO  EL 

As  a  aeans  of  coapariscn  of  aerosol  initialization 
schemes,  the  sodel  cffars  the  option  of  initialization  with 
equilibrium  values.  These  initial  values  are  based  on  data 
collected  fzos  the  J1SXN  Experiment  {Fairall,  et  al,  1982}  . 
A  large  aercsol  data  base  was  collected  during  equilibrium 
ccnditicrs,  with  respect  to  wind  speed,  during  12  hour 
periods  and  then  normalized  to  a  reference  relative 
humidity.  Aerosol  spectra  were  then  grouped  into  six 
different  wind  speed  ranges.  Fcr  the  equilibrium  initialize' 
tion  scheme  in  the  NFS  sodel,  a  pre-assigned  value  cf  the 
five  radii  are  used  as  initial  aerosol  values,  based  on  the 
wind  speed  at  the  model  start  time.  An  example  of  the  equi- 
libriui  aercsol  spectra  from  the  JASIN  Experiment  is  shewn 
in  figure  2.4.  mote  that  the  reference  relative  humidity  for 
this  experiment  was  87  percent.  The  graph  is  referenced  here 
only  tc  show  an  equilibrium  type  distribution;  the  initial 
equilibrium  and  HAGAT  aerosol  values  will  be  shown  in 
chapter  IY.  Once  initialized  with  either  equilibrium  data  or 
calculated  data,  the  model  will  run  in  the  same  manner. 


figure  2-4  Aercsol  spectra  as  a  function  of  size 
for  selected  wind  speeds. 


This  equilibrium  scheme  illustrates  the  fundamental 
difference  between  the  NFS  model,  which  is  a  dynamic  model, 
and  a  steady-state  type  model.  A  representative  mcdel  using 
this  equilibrium  approach  is  the  WMK  model,  mentioned  in  the 
beginning  of  this  chapter.  Ihis  type  of  output  describes  the 
marine  aerosol  distributions  as  a  function  of  wind  speed 
(surface  generation)  ,  relative  humidity  (growth  factor)  ,  and 
elevation  (vertical  variation  of  aerosols  with  height 
assuming  a  steady  state  vertical  transport  process) .  The 
mathematical  representation  of  this  model  /Fairall,  et  al, 
1982a ]  is 

n(r)  =  (r/a)*1.62  (c1  +  c2vi)F*exp  (-z/,h0F-8.5(r/a)  ^  (2.5) 


where 

r*  the  particle  size  in  microns, 

u*  tie  wind  speed, 

v*  0.5  for  u  S  4  m/s, 

v*  u  -3.5  for  u  >  4  m/s, 

f*  U(v/60)3, 

r*  o. 3 ea-o. 0029  3V1”  , 

Z*  height  above  sea  surface,  m  , 
h  »  scale  height, m  (800  m  for  Z<1000  m)  , 
a*  O.dlezp  (0.  066s/ (1. 0 58-s) )  , 

S*  H/ 100,  (H*  relative  humidity  in  percent). 

The  ether  constants  are  given,  based  on  the  value  of  v. 


▼,  «/s 
v  <  7 
v  >  7 


1000 

6900 


1.15 
0.  29 


The  equation  described  produces  a  number  density  spectrum, 
n(r)  .  This  relates  to  the  volume  spectrum  V(r),  as  shewn  in 
equation  2.2. 


III.  £AJA  AC0PISITI08  ABD  SYNOPTIC  CONDITIONS 
1.  DATA  ACCOISITIOH 

The  data  were  obtained  from  a  ship  and  aircraft  experi¬ 
ment,  flaiiflf  Aerosol  generation  and  Transport  (MAGAT),  held 
from  28  April  to  9  Say  1980,  in  a  region  30  to  50  miles  cff 
the  coast  cf  Monterey  Bay,  Ca.  Data  gathering  was  done  with 
instruments  mounted  cn  both  an  aircraft  and  a  surface  ship, 
the  B/V  ACANIA.  Aerosol  data  were  obtained  from  the  Airborne 
fiesearch  Associates  turbo-charged  Bellanca,  using  a  Farticle 
Measuring  System  (PHS)  Axial  Symmetric  Scattering  Aerosol 
Probe  (ASSAP)  particle  counter.  All  measured  data  were 
sampled  every  2.5  seconds,  with  a  two-scan  average  cf  every 
five  seconds.  The  scans  were  collected  in  "ladder"  profiles, 
during  which  the  aircraft  made  measurements  at  a  constant 
altitude  fcr  two  minutes,  climbed  to  a  new  altitude,  and 
repeated  a  rew  measurement  run.  The  instrument  utilized  60 
size  channels  from  2.8  to  14.0  micron  radius.  In  most  cases 
the  ladders  extended  from  near  the  sea  surface  (3  m)  up  up 
through  the  well-mixed  boundary  layer,  to  a  few  steps  above 
the  inversion.  The  elevation  was  generally  up  to  5  kilome¬ 
ters.  A  typical  ladder  profile  contains  10  to  14  steps.  The 
step  heights  were  randomly  chosen,  but  an  attempt  was  made 
to  keep  each  step  height  consistent  between  ladders.  Air  and 
dew-point  temperatures  were  also  measured  and  used  to  calcu¬ 
late  relative  humidity  for  the  correction  factor  mentioned 
in  chapter  two.  The  aircraft  also  flew  ascending  spirals  (in 
the  vicinity  of  the  ladders)  during  which  other  meteorolo¬ 
gical  parameters  were  collected.  This  data  yielded  vertical 
soundings  similar  to  those  provided  by  radiosondes. 


There  were  two  aerosol  instruments  fitted  on  the  R/V 
ACANIA.  One  was  the  PMS  model  CSAS  (classical  scattering) 
and  the  ether  was  model  ASAS  (active  scattering) ,  controlled 
ty  a  INS  data  acquisition  systaa  (DAS-32)  with  a  computer 
interface.  The  shipboard  system  measured  aerosols  in  90 
different  size  channels  from  0.9  to  14.0  micron  radius. 

It  was  ncted,  during  the  flybys  over  the  R/V  ACANIA,  that 
the  aircraft  measurements  did  not  agree  with  that  of  the 
ship.  The  size  distributions  from  the  aircraft  measurements 
were  consistently  smaller  than  the  ship  measurements,  for 
radii  greater  than  1.C  micrcn.  In  addition,  the  differences 
increased  with  radius.  Since  the  shipboard  aerosol  system  is 
newer  and  had  a  wider  size  range,  and  better  sensitivity, 
aircraft  aerosol  data  were  corrected  to  agree  with  the  ship 
aerosol  data  /Fairall,  1980  and  Fairall  et  al,  198_q7  . 

Profiles  of  virtual  potential  temperature  and  mixing 
ratic  were  obtained  from  three  different  sources.  The  reason 
for  using  these  paraceters  instead  of  temperature  and  dew 
point,  is  that  mixed  layer  inversions  are  more  easily  iden¬ 
tified  with  the  former  variables.  The  sources  were  the 
spiral  flights  from  the  aircraft,  the  radiosonde  launches 
from  BPS  and  from  the  R/V  ACANIA. 


B.  SIBOETIC  DATA 

Surface  and  500  ib  synoptic  charts  and  the  GOES  1EST 
satellite  images  were  used  to  evaluate  the  synoptic  condi¬ 
tions.  Charts  are  from  the  NOAA  weekly  series  cf  daily 
weather  maps.  In  addition  to  data  collected  from  the  R/V 
ACANIA  and  the  aircraft,  local  weather  data  were  also  avail¬ 
able  fxca  the  0. S.  Arty's  Pritzsche  Field  weather  facility. 
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ligura  3.2  SOES  R«st  satellite  iaagery,  0915  POT  2  Hay,  1980. 


Figure  3*4  Saae  a a  Figure  3,2  except  1645  PDT ,  a  Bay 


Figure  3.6  Same  as  figure  3.2  except  1745  PDF,  5  Bay. 
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Figure  3.1  Saee  as  figure  3.2  except  1645  PDT ,  6  Bay. 


lort  Crd,  Ca.  Satellite  photographs  and  NO AA  aaps  for  the 
two  different  tine  periods  <2-4  Bay  and  6-7  May)  ars  shewn 
in  figures  3.1  through  3.8.  Each  three-  day  sequence  depicts 
the  24  hour  period  being  considered;  24  hours  before  and 
after.  Atmospheric  soundings  for  the  two  model  runs  are 
shown  in  chapter  IV,  figures  4.1  and  4.2. 

C.  5IN0ETIC  CONDITIONS 

Several  weak  frontal  systems  passed  through  th9  area 
during  the  experiment.  Showers  occurred  during  the  first  and 
last  cays,  associated  with  the  fronts.  Low  cloudiness  and 
fog  cccured  during  the  morning  from  29  April  to  5  Hay,  with 
fog  returning  again  cn  the  last  day. 

At  the  beginning  cf  the  period,  the  area  was  dominated  by 
a  slcwly  eastward  migration  of  a  cut-off  low  at  500  mb.  By 
early  morning  on  2  (lay,  the  area  was  under  the  influence  of 
a  weak  ridge.  On  3  Bay  the  area  was  under  divergent  flew  at 
the  upper  level.  An  upper  level  low  had  formed  off  the  Baja 
Ca.  coast  cn  4  Bay,  leaving  the  area  under  the  influence  of 
a  col.  Cn  5  Hay  the  area  was  between  an  upper- level  trough 
and  ridge,  and  by  6  Bay  the  area  was  on  the  back  side  of  the 
trough.  Eecause  cf  the  deepening  of  the  trough,  the  area  was 
still  cn  the  back  side  of  the  trough  on  7  Hay.  A  new  upper 
level  treugh  formed  and  approached  the  area  on  the  final  two 
days  cf  the  experinert. 

Surface  winds  were  relatively  light,  0  to  10  kts,  at  the 
beginning  of  the  peried,  and  increased  toward  the  end  of  the 
peried  tc  16  kts,  gusting  to  22  kts. 

An  ixportant  feature  of  these  interpretations  is  the 
nature  of  the  mixed-layer,  often  topped  by  an  inversion, 
with  regcard  to  stability,  and  therefore  mixing  intensities. 
It  is  assumed  that  the  mixing  becomes  greater  as  conditions 
become  acre  unstable. 
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It  tke  start  cf  the  expedient,  the  nixed  layer  exhibited 
stable  tc  slightly  unstable  conditions,  until  about  1800 
local  tine  (POT)  28  April,  when  conditions  became  acre 
neutral.  The  neutral  condition  remains  until  1  Hay  when 
conditiccs  cnce  again  become  stable.  A  weak  frontal  passage 
before  0500  PDT  on  25  April  does  not  appear  tc  affect  the 
■ixed  layer  profile.  The  layer  remains  stable  until  a 
frontal  passage  on  2  Hay  when  conditions  beccae  neutral  and 
resales  sc  cntil  5  Hay.  On  the  lorning  of  5  Hay,  conditions 
are  slightly  stable,  but  returned  to  neutral  on  6  Hay, 
despite  a  frontal  passage  at  1300  POT  on  5  Hay,  and 
regaining  neutral  to  the  end  of  the  experiment  on  9  Hay. 

0.  Sill  SI1BCTICI 

Data  free  the  HAGAT  experiment  were  selected  for  model 
verification  on  the  basis  of  the  relative  positions  of  the 
two  data  gathering  platforms  and  the  general  synoptic  ccndi- 
tions.  The  locations  of  the  surface  ship  and  the  aircraft 
did  not  always  coincide  during  the  experiment.  Therefore 
soundings  and  aerosol  data  were  considered  suitable  if  the 
two  were  within  20  nautical  miles.  Hourly  wind  speed, 
direction,  and  SST  were  gathered  from  the  ACANIA  and  were 
used  along  with  the  aircraft  soundings.  Plots  of  wind 
speed,  air  (solid  line)  and  sea-surface  (dashed  line)  temp¬ 
erature,  and  relative  humidity  for  the  two  24  hour  periods 
are  shown  in  figures  3.9  and  3.10.  Another  consideration  of 
the  aircraft  data  was  the  relative  locations  of  the  aerosol 
ladders  and  the  soundings.  Although  both  gathered  during 
vertical  flight  profiles,  they  were  not  done  at  the  same 
time.  The  synoptic  conditions  during  the  experiment  include 
perieds  cf  frontal  weather  which  was  also  avoided  in  the 
analysis.  Two  time  blocks 


RELATIVE  HUMIDITY  C4)  TEMPERATURE  (C)  WIND  SPEED  (M/S) 


FLT#6 
3  MAY  80 


Figare  3.11  Flight  p^tl^of  ajjcraft  on  the  aorning  of 
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FLT#8 
4  MAY  80 
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FLT  #14 
7  MAY  80 


Hgnre  3.16  sill  it  figure  3.13  »ei;t  afternoon  of  7  Hay  I960 


Ilgara  3.17  Crcisa  track  of  B/7  &C&NI&  3  Ha 7  1980. 
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3.16  except  4  Hay  1980 


were  chosen;  3-4  May  and  6-7  May.  Figures  3.11  through  3.16 
indicate  the  routing  of  the  aircraft  flights  for  the  two 
tiae  blocks,  and  figures  3.17  through  3.19  depict  the  course 
cf  the  B/V  AC A N I A  daring  the  sane  times. 

Atmospheric  data  from  the  aircraft  soundings  were  used 
instead  cf  the  radicscndes  from  the  R/V  ACA NIA .  This  was 
done  for  two  reasons:  first  because  the  aircraft  soundings 
contained  acre  levels,  and  second  they  ware  in  closer  prox¬ 
imity  tc  the  aircraft  aerosol  ladders.  Atmospheric  and 
aerosol  data  for  the  first  run  was  taken  from  flight  6,  on 
the  acrning  of  3  May.  The  locations  and  times  are  shown  in 
figure  3.9.  The  aircraft  was  approximately  45  miles  from  the 
coast  when  a  sounding  was  taken  at  1143  PDT,  and  an  aerosol 
ladder  (111)  was  taken  at  1200  PDT.  Data  for  the  two  verifi¬ 
cations  of  atmospheric  parameters  and  aerosol  values  were 
chosen  based  on  their  proximity  to  the  initial  data.  A 
sounding  and  aerosol  ladder  on  the  afternoon  of  3  May 
(sounding  time  of  1714  PDT,  and  aerosol  ladder  113)  approxi¬ 
mately  eight  hours  frcm  the  initial  time  was  chosen  for  the 
first  verification,  and  is  depicted  in  the  route  of  flight 
7,  figure  3.10.  Data  for  the  second  verification  was  chosen 
from  flight  8,  figure  3.11;  sounding  time  1151  PDT,  and 
aeroscl  ladder  116.  Data  for  the  two  verifications  are  bcth 
within  20  miles  of  the  initial  data. 


Data  for  the  second  model  run  is  shown  in  figures  3-12 
through  3-14.  The  model  was  initialized  with  data  from  the 
>732  PDT  scunding  and  ladder  122,  approximately  25  miles  eff 
the  Hcnterey  coast.  The  first  verification  data  was  taken 
from  flight  12  the  next  morning  approximately  17  hours 
later;  sounding  tiae  1743  PDTand  ladder  123.  The  seccnd 
verificaticn  data  was  taken  from  flight  13,  approximately  24 
hours  frcm  the  initial  data;  sounding  time  1836  PDT  and 
ladder  124.  Data  collected  for  the  two  verifications  were 
within  5  tiles  of  the  location  of  the  initial  data. 


Additional  wind  speed,  direction,  and  SST  collects  from 
the  B/V  AC ANIA  were  used  for  initialization  and  verification 
cf  the  two  model  runs.  As  mentioned  in  the  model  descrip¬ 
tion,  up  to  10  forecast  values  for  each  can  be  incut  into 
the  model.  However,  due  to  the  continuous  movement  of  the 
B/V  ACANIA  during  the  MAGAT  experiment,  where  SST  values 
were  charging  because  of  strong  coastal  gradients,  a  single 
value  of  SST  was  used.  This  value  was  chosen  when  the  fi/V 
ACANIA  was  closest  to  the  areas  where  the  aircraft  data  was 
taken  fcx  each  cf  the  two  initialization  sites.  The  wind 
speed  and  direction  cata  from  the  B/V  ACANIA  was  considered 
representative  of  tHe  data  collection  sites,  and  therefore 
input  into  the  model  runs  at  approximately  3  hour  intervals. 


IT.  HOD El  IBITI ALIZATION  AND  BESDLTS 


1.  DIGITIZED  ATHOSPEEBIC  SOOBEIHGS 

The  atacspheric  inputs  for  the  model  (chapter  2)  are: 

1.  The  nixed  layer  equivalent  potential  temperature  (C) 
and  specific  humidity  (gm/kg). 

2.  The  "jump"  disccrtinuities  in  each  of  the  above  values. 
The  junp  meaning  the  difference  between  the  mixed  layer 
value  and  the  value  at  the  top  of  the  inversion. 

3.  The  lapse  rate  for  each  parameter  value  atcve  the 
inversion. 

4.  The  depth  of  the  well  mixed-layer. 

To  simplify  the  process  of  calculating  these  values  from  the 
soundings,  a  digitizing  scheme  was  designed  by  the 
Environmental  Physics  Group  at  the  Naval  Postgraduate 
School.  This  scheme  transforms  the  sounding  into  the  struc¬ 
ture  that  describes  the  necessary  model  inputs.  Examples  are 
shown  in  figures  4.1  and  4.2.  Note  that  the  "jumps"  are 
depicted  as  occuring  in  an  infinitely  thin  layer. 

E.  ADDITIONAL  HOD EL  INPOT  AID  ADJOSTHENTS  TO  THE  DIGITIZED 
SOUNDINGS 

The  reascn  for  making  adjustments  to  the  atmospheric  part 
cf  the  model  is  to  provide  a  better  basis  for  evaluating  the 
behavior  of  the  aeroscl  prediction.  This  is  optimized  when 
the  mcdel  is  producing  the  best  possible  forecast  cf  the 
atmospheric  parameters.  The  parameters  are  forecast  very 
well  in  both  model  runs,  with  these  adjustments  made.  The 
model  is  accurately  describing  the  height  of  the  inversion 
through  adjustments  in  the  subsidence  rate,  and  the  lifting 
condensation  level  (LCL)  is  accurately  generating  clouds 


through  adjustments  tc  the  temperature  and  specific  humidity 
which  are  within  the  limits  of  instrumentation  error.  The 
cloud  verifications  were  made  with  satellite  data  and  obser¬ 
vations  from  the  R/V  1CANIA  log. 

An  initial  input  for  the  model  is  the  subsidence  rate. 
Using  a  first  guess  cf  -.005  m/sec,  the  model  generated  a 
plot  cf  the  inversion  height  for  the  24-hour  period.  Eased 
on  the  verification  data,  the  subsidence  rate  was  further 
adjusted  tc  bring  both  the  predicted  and  observed  inversion- 
height  values  together.  Then  a  small  adjustment  was  made  to 
the  mixed  layer  values  to  match  the  stability  and  cloud/ 
cloud  free  patterns  with  what  was  actually  observed  during 
the  period.  In  the  first  model  run  (3  May)  ,  the  potential 
temperatcre  was  increased  0.5  degrees  and  the  specific 
humidity  was  decreased  1  gm/kg.  The  resultant  subsidence 
rate  necessary  tc  match  the  predicted  and  observed  inversion 
height  was  -.0042  m/sec.  These  adjustments  yield  an  atmos¬ 
pheric  profile  similar  to  what  was  actually  observed.  The 
air  mass  became  slowly  saturated  enough  to  generate  clouds 
in  the  upper  part  of  the  mixed  layer  at  about  0300  PDT  on  4 
May.  Piets  of  the  irversion  and  LCL  behavior,  and  the  temp¬ 
erature  and  humidity  profiles  are  shown  in  figure  4.3.  In 
the  plot  fer  LCL  and  inversion  height,  the  letter  WHH 
depicts  the  verification  times  and  values  for  the  inversion 
height,  and  the  letter  wLn  similarly  indicates  values  for 
the  I CL  verifications.  The  letters  "I"  and  "QM  in  the  temp¬ 
erature  and  specific  humidity  graphs  indicate  verification 
values  fer  each  respective  value. 

The  atmcspheric  plots  for  the  second  model  run  (6  May) 
are  seen  in  figure  4.4.  For  this  run  the  specific  humidity 
was  reduced  from  6.4€  to  6.3  gm/kg,  the  potential  tempera¬ 
ture  was  unchanged,  and  a  subsidence  rate  of  -0.01  m/sec  was 
input  for  the  run. 


Ilguts  4.1  Digitized  sounding  for  first  sodel  run, 
3  Bay  1980,  1143  PDT. 


49 


RECORD  ♦  33 

/ 


W9«.  ».2  *«vk?ifliS.<,i?<9*fg;’t  £ot  s,eo»4 1 


Height  (km)  0  (C) 


lo  | — r~r  t  -i  t— j  r  r  t  r  t-{  t  r  r*r  r  j  r  t  r~r—r 


O) 

\ 

CD  8 


18 


T  i  me 


18 


Fignra  4.4  Sana  as  Figure  4.3,  azcapt  for  second  run 


C.  A1B0SC1  IHITIALI2STIOH  EBOCEDURES 

The  aerosol  data  were  gathered  in  ladders,  usually  up  to 
five  kilometers,  and  were  of  actual  concentrations  of  parti¬ 
cles  for  a  given  increment  of  radius,  n  (r)  .  This  number 
contained  bcth  marine  and  continental  components,  collected 
at  the  axbient  relative  humidity. 

The  initialization  scheme  was  performed  in  three  main 
steps.  The  first  step  was  to  identify  the  marine  and  conti¬ 
nental  components.  By  using  the  digitized  soundings  for  each 
model  run,  a  determination  was  made  of  two  height  values: 
one  representing  the  well  mixed-layer  and  the  other  repre¬ 
senting  the  free  troposphere,  above  the  inversion.  Both  the 
mixed  layer  and  free  troposphere  aerosol  counts  for  each  of 
the  five  radii  were  calculated.  The  second  step  was  to  apply 
the  humidity  growth  factor  to  the  values  and  calculate 
dv/dr.  The  relationship  between  n  (r) ,  which  was  measured, 
and  dv/dr  is  given  by  equation  2-2.  The  final  step  was  to 
subtract,  fcr  each  radius,  the  free  troposphere  dv/dr  value 
from  the  nixed  layer  dv/dr  value.  fiith  the  assumption  of  a 
constant  prefile  of  continental  aerosols  through  both  the 
nixed  layer  and  into  the  free  troposphere,  this  subtraction 
leaves  in  only  the  mixed  layer  marine  components,  normalized 
to  80  percent  relative  humidity.  The  values  of  dv/dr  used 
to  initialize  the  mcdel,  both  equilibrium  and  HAGAI,  are 
shown  in  tables  4.1  and  4.2. 

This  three  step  prccess  was  dene  for  both  model  runs,  and 
for  two  verification  times  within  each  run.  As  described  in 
chapter  III,  the  verification  ladders  were  chosen  to  be  as 
close  as  possible  tc  the  position  of  the  initial  aerosol 
data.  In  both  cases,  the  verification  data  were  within  15 
miles  of  the  initial  data.  The  plots  of  the  aerosol  output. 


toth  for  the  HAGAT  and  equilibrium  initializations,  is 
depicted  as  log  10  (dv/dr)  against  the  24  hour  time  period, 
in  half  hour  time  steps.  In  both  (initialization)  cases,  the 
WX"  values  represent  the  HAGAT  verification  values.  In  both 
model  runs,  aerosol  output  for  the  HAGAT  initialization 
schese  is  represented  by  a  continuous  line,  and  the  dashed 
line  depicts  the  run  for  the  equilibrium  initialized  cases. 

T1BLE  I 

Initial  equilibrium  aerosol  values 
radius,  urn 


o.e 

2.0 

5.0 

10 

V 

6.  C 

2.8 

3.  4 

3.0 

0.3 

i 

9.0 

10 

10 

5.0 

1.0 

n 

11 

18 

25 

13 

6.4 

a 

13 

20 

30 

15 

15 

15 

33 

22 

22 

22 

s/s 

18 

28 

35 

28 

28 

TABLE  II 

Initial  observed  aerosol  values 
radius,  ua 


0.  E 

2.0 

5.0 

10 

3  Hay  run 

42.2 

21.8 

2.8 

0 

6  Hay  run 

45.4 

22.8 

5.2 

5.1 

Whet  considering  the  equilibrium  initial  values  from 
table  one,  it  is  easy  tc  see  the  relationship  between  wind 
speed  and  production.  In  almost  all  cases,  the  aerosol 
volume  counts  increase  as  the  wind  speed  increases.  The 
distribution  of  radii  for  a  given  wind  speed,  however,  is 
not  as  well  defined.  For  lower  wind  speeds,  the  larger 
radii  are  few  it  nuiber,  and  as  the  wind  speed  increases 
there  is  a  dramatic  increase  in  the  number  of  larger  parti¬ 
cles.  This  is  reasonable  due  to  white  cap  production.  Also, 
there  appears  to  be  a  slight  maximum  of  particles  at  2.0 
microns,  reguardless  of  wind  speed.  The  15  micron  initiali¬ 
zation  values  in  the  tables  and  the  model  output  fcr  the 
same  size  are  not  shown.  This  is  because  the  observed 
initial  values  collected  from  the  aircraft  were  zero. 

In  comparing  the  EAGAT  initial  values,  one  finds  that  the 
smaller  radii  are  approximately  twice  the  equilibrium  count 
for  even  the  largest  wind  speed.  It  is  interesting  tc  note 
that  for  such  a  large  variation  in  initial  wind  speed  for 
the  two  days,  there  is  very  little  variation  in  the  initial 
values,  for  all  the  the  particle  sizes.  Based  on  results  of 
Bonahan  /et  al  1983/  ,  it  was  also  surprising  to  find  no 
volume  counts  for  the  larger  particle  sizes  since  wind  speed 
was  above  10  a/sec.  For  the  strong  wind  case,  there  were  no 
aerosols  above  10  microns  collected  from  the  aircraft. 

B.  I IB ST  BCDEL  SOB  F1SULTS 

The  plots  of  dv/dr  over  the  24  hour  period  for  the  first 
model  run  are  shown  in  figures  4.5  and  4.6.  Fcr  each  radius, 
the  solid  line  indicates  the  profile  with  the  NAGAT  initial¬ 
ized  values,  and  the  dashed  line  represents  the  profile  with 
the  equilibrium  initialization. 


the  model  shows 


For  the  three  smaller  sizes. 


change  ever  the  period,  however  the  model  does  produce  a 
trend  in  aerosol  production.  ahatever  volumes  were  used  to 
initialize,  the  number  tends  to  stay  within  one  order  of 


magnitude  ever  the  period.  The  verification  values  compare 
very  well  to  the  model  run  with  the  HAG AT  initialization. 
The  model  and  verification  dv/dr  values  differ  by  less  than 
10  units.  For  the  larger  size,  10  microns,  note  the  quick 
production  of  aerosol  when  initialized  with  with  a  zero 
value.  Bithin  four  hours,  the  model  is  producing  close  to 
the  same  values  as  the  equilibrium  case.  From  this  feur-heur 


point,  the  two  cases  behave  very  much  the  same.  There  is  a 
diurnal  variation  evident  in  the  larger  size  cases;  a  lag  is 
noticeable  in  the  prcduction  through  the  night  time,  with  as 
increase  after  sunrise  and  through  the  morning.  The  verifi¬ 
cation  values  for  these  larger  sizes  are  not  good  for  either 
verification  scheme. 
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Pignre  4.6  Saie  as  fig ora  4.5,  except  a)  r«  5 

and  6)  r*  10  an. 


E.  S  ICC  HD  HODEL  HOE  BESOLTS 


Aeroscl  plots  for  loth  the  equilibrium  and  MAG AT  initial¬ 
ization  schemes  are  shown  in  figure  4.7  and  4.8.  The  dv/dr 
values  for  both  schemes  are  higher,  due  mainly  to  the 
increased  wind  speed.  For  all  of  the  particle  sizes,  there 
is  mere  cf  a  diurnal  variation.  This  is  slightly  evident  in 
the  smaller  radii,  and  veil  defined  in  the  larger  radii. 
This  variation  is  strengthened  by  a  wind  speed  minimum  just 
before  sunrise.  Note  also  that  the  drop  in  production  of 
all  the  radii  around  0000  PCT  coincides  with  both  the  slight 
drop  in  wind  speed  from  10  tc  8  m/sec,  and  the  formation  of 
clouds,  evident  in  the  LCL  and  inversion  figure.  This 
decrease  in  production  due  tc  cloud  formation  is  not  notice¬ 
able  in  the  first  model  run.  The  two  initialization  schemes 
are  ic  tore  agreement  than  with  the  first  run,  again  due 
mainly  tc  the  higher  wind.  With  the  larger  wind  speeds  in 
this  run,  the  model  is  generating  values  closer  to  the  veri¬ 
fication  values  (in  teth  initialization  schemes)  than  in  the 
first  rue,  where  the  wind  speed  is  significantly  less. 


I.  SUHH1BI  OF  RESULTS 

The  model  is  producing,  with  a  minimus  of  adjustment,  an 
accurate  description  of  the  atmospheric  boundary  layer 
parameters.  Once  a  subsidence  rate  is  selected  such  that  the 
24-hcur  predicted  inversion  height  approaches  the  observed 
value,  the  values  of  temperature  and  humidity  only  require  a 
minimum  cf  adjustment.  In  both  cases  the  verification  values 
cf  temperature  were  within  2  *  C,  and  the  values  of  specific 
humidity  within  1  gm/kg.  More  significantly,  the  trends 
ever  the  24-hour  periods  for  both  temperature  and  humidity 
were  accurately  described  by  the  model,  as  defined  by  the 
verification  data. 

With  respect  to  the  aerosol  input,  the  model  is  also 
producing  the  correct  trend  in  behavior.  There  is  evidence 
in  the  output  for  the  trend  in  aerosol  production  to  be  a 
function  of  not  only  wind  speed  and  relative  humidity,  but 
also  the  concentration  of  the  initial  aerosol  values.  This 
is  seen  in  the  first  run,  particle  sizes  0.8  and  2.0 
micrcrs.  The  atmospheric  parameters  are  the  same  in  both 
cases,  for  each  radius;  but  it  is  the  initialization  scheme 
that  generates  a  different  trend  in  each  case.  The  model  run 
with  the  HA  GAT  initial  values  produces  a  loss  in  production, 
and  the  verification  values  are  within  one  order  of  magni¬ 
tude  cf  the  model  values,  with  the  equilibrium  case,  the 
trend  is  an  increase  in  production  over  the  time  period, 
with  the  model  verification  values  two  orders  of  magnitude 
away  from  the  observed  verification  values. 

The  production  trends  in  the  second  model  run,  with  both 
initialization  schemes,  are  much  the  same.  In  tk\s  case 
there  is  less  of  a  difference  between  the  model  verification 
values  fer  the  two  initializations.  with  the  exception  of 
the  smallest  radius,  0.8  microns,  the  verification  values 


for  all  tfae  radii  of  the  second  run  are  vithin  one  order  of 
sagnitude  cf  the  model  output.  The  fact  that  the  MAGAT 
initial  values  are  closer  tc  an  equilibrium  value  at  a 
higher  wind  speed  seems  tc  imply  a  less  sensitive  model, 
when  initialized  with  a  higher  wind  speed.  Accordingly, 
eeroscl  production  is  high;  and  the  model  accurately 
predicts  the  behavicf  of  the  aerosol  in  the  stronger  wind 


V.  COHCIOSIONS 


The  boundary  layer  model,  with  a  minimum  cf  input,  will 
produce  an  accurate  description  of  the  MABL  [Davidson  et  al, 
1984 ]  and  the  aerosol  behavior  within  the  marine  boundary 
layer.  The  model  is  generating  values  for  marine  particles 
within  this  layer,  and  whenever  there  is  a  large  difference 
between  the  predicted  and  observed  values  ,  it  was  the 
observed  values  that  were  toe  high.  This  observation 
suggests  that  the  observed  values  collected  by  the  aircraft 
(made  up  of  marine  and  continental  particles)  may  have 
included  mere  continental  than  originally  expected.  A 
linear  distribution  (vertically,  through  the  boundary  layer) 
was  assumed  in  the  aerosol  initialization  process. 

Here  significantly,  the  model  is  correctly  identifying 
the  production  trend  cf  aeroscls  over  the  time  period.  The 
model  is  run  with  identical  atmospheric  inputs,  and 
different  aerosol  initialization  schemes  produce  two 
different  production  trends.  Tne  wind  speed  and  relative 
humidity  (the  only  atmospheric  inputs  of  previous  models) 
are  the  sane,  and  the  model  is  generating  different  profiles 
based  on  initial  aercscl  concentrations.  Bith  a  minimum  of 
two  case  studies,  it  is  clear  that  the  process  of  entrain¬ 
ment  cf  air  from  above  the  mixed  layer,  containing  nc  marine 
particles,  is  affecting  the  behavior  of  the  model  output. 
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